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We find that scale-free random networks are excellently modeled by a deterministic graph. This 
graph has a discrete degree distribution (degree is the number of connections of a vertex) which is 
characterized by a power-law with exponent 7 = 1 + In 3/ In 2. Properties of this simple structure 
are surprisingly close to those of growing random scale-free networks with 7 in the most interesting 
region, between 2 and 3. We succeed to find exactly and numerically with high precision all main 
characteristics of the graph. In particular, we obtain the exact shortest-path-length distribution. 
For the large network (IniV 3> 1) the distribution tends to a Gaussian of width ~ Vhi N centered 
at t ~ IniV. We show that the eigenvalue spectrum of the adjacency matrix of the graph has a 
power- law tail with exponent 2 + 7. 
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The essence of the modern situation in network sci- 
ence is the change-over from study of classical random 
graphs with Poisson degree distributions |I]J^] to explo- 
ration of complex networks with fat-tailed degree distri- 
butions |?]-^) . The prominent particular case of such nets 
are networks with power-law degree distributions (scale- 
free networks) ||. While growing, such nets actually self- 
organize into scale-free structures. These networks play 
a great role in Nature The Internet, the WWW, 

and many basic biological networks belong to this class. 
Fat tails of the degree distributions produce a number of 
intriguing effects jg-[L4j. 

Such networks are widespread, but very little is still 
known even about their basic properties 15 lq|. Most 



of real growing scale-free networks have 7 exponent of 
the degree distribution P(k) ~ k 1 in the range (2,3), 
but this case turned to be the most difficult and unex- 
plored. In particular, no exact results for the average 
shortest-path length I are known in this situation. The 
only known exact shortest-path length distributions were 
obtained for the simplest equilibrium networks [|l7l . [l"8i . 
Notice that if 7 < 3, standard estimates of J |19| are 
inapplicable to equilibrium networks with uncorrelated 
vertices. Correlations in growing networks are inevitable, 
and the results are even less encouraging. The generic 
property of these networks, which makes their analytical 
study so hard, is a complex structure of their adjacency 
matrices. 

Scale-free random networks naturally have a contin- 
uous degree distribution spectrum, but it has recently 
been shown that discrete degree distributions of some de- 
terministic graphs also have a power-law decay ]2(J . The 
leading idea of the present work is the following. If it is so 
hard to get exact or precise results for scale-free random 
networks, especially, with 2 < 7 < 3, let us (i) construct 
a simple deterministic scale- free graph with such 7, (ii) 
obtain exact (analytical) and precise (numerical) answers 
for main structural and topological characteristics of the 
graph, (iii) compare its properties with known charac- 



teristics of corresponding scale-free growing random net- 
works, and (iv) if these properties coincide, make new 
predictions for scale-free growing random networks. We 
actually model scale-free networks with 2 < 7 < 3 by the 
deterministic graph, whose structure can be described 
completely. 

Here we present results of this program. We succeed 
to find a number of exact characteristics of the scale-free 
deterministic graph, which are still unknown for random 
scale-free networks. The structural properties of deter- 
ministic and random scale-free growing networks proved 
to be surprisingly close to each other, so that our results 
can be reasonably applied to random growing nets. 

Pseudofractal graph. — The growth starts from a single 
edge connecting two vertices at t = — 1 (see Fig. |l|). At 
each time step, to every edge of the graph, a new vertex 
is added, which is attached to both the end vertices of 
the edge. Then, at t = 0, we have a triangle of edges 
connecting a triple of vertices, at t — 1, the graph con- 
sists of 6 vertices connected by 9 edges, and so on. The 
total number of vertices at "time" t is N t = 3(3* + l)/2, 
and the total number of edges is L t — 3 t+1 , so that the 
average degree is k t = 2L t /N t — 4/(1 + 3 - '). 

This simple rule produces a complex growing network 
which is certainly not a fractal plfl . Indeed, at any step, 
the entire graph can be set inside of a unit triangle. 
This means that the structure has no any fixed finite 
fractal dimension. On the other hand, one can depict 
the graph in another way, namely, like in Fig. [l] where 
the graph is surrounded by a long chain of edges. The 
length of this "perimeter" is Pt = 3-2* edges, whence 
N t ~ L t ~ p' n3 / ln2 . We can, however, fold the graph 
into a more compact structure with a different border, 
so that In 3/ In 2 is not a fractal dimension of the struc- 
ture but only some characteristic value. We failed to 
introduce a well defined spectrum of fractal dimensions 
(between In 3/ In 2 and 00), hence the network cannot be 
called a multifractal. Thus, this graph is not a fractal 
but only parody of it, and we call it, for brevity, pseud- 
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ofractal. Notice that the graph contains numerous loops 
and hence is very far from tree-like. 

Adjacency matrix. — By definition, an element of an 
adjacency matrix is equal to 1 or depending on whether 
an edge between vertices i and j is present or not. The 
adjacency matrix A t structure is schematically shown in 
Fig. 0. At t = — 1, this is the 2x2 matrix with zeros 
on the diagonal and two unit elements. At timestep i, 
we add rows and columns i, j = AT t _i + 1, . . . , Nt (new 
vertices) to the matrix. Matrix is symmetric, and each 
unit element above the diagonal of the matrix A t -\ 
generates, in addition, two unit elements ai S and aj S of 
At- Here Nt-i + 1 < s < Nt. Other elements are zeros. 
This produces the sparse block matrix shown in Fig. 

Degree distribution. — The degree spectrum of the 
graph is discrete: at time t, the number m(k, t) of ver- 
tices of degree k = 2, 2 2 , 2 3 , . . . , 2*" 1 , 2*, 2 t+1 is equal to 
3', 3 t_1 , 3 t_2 , . . . , 3 2 , 3, 3, respectively. Other values of 
degree are absent in the spectrum. Clearly, for the large 
network, m(k,t) decreases as a power of k, so the net- 
work can be called "scale-free" . Spaces between degrees 
of the spectrum grow with increasing k. Therefore, to 
relate the exponent of this discrete degree distribution to 
standard 7 exponent of a continuous degree distribution 
for random scale- free networks, we use a cumulative dis- 
tribution Pcum(k) = Ylk'>k m (k', t)/Nt ~ fc 1 ^ 7 . Here k 
and k' are points of the discrete degree spectrum. Thus 
we obtain 



In 3 
^ = 1+ hT2' 



(1) 



so that 2 < 7 = 2.585 ... < 3. Compare 7 with 
the characteristic exponent in the relation between the 
"mass" and the "perimeter" of the graph. Also, no- 
tice that the maximal degree of a vertex is equal to 
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which coincides with a 



standard relation for the cutoff of degree distribution in 
growing scale- free networks || . 

Distribution of clustering. — By definition, the cluster 
coefficient C of a vertex is the ratio of the total number 
of existing connections between all k its nearest neigh- 
bors and the number k(k — l)/2 of all possible connec- 
tions between them. Usually, only the average value of 
the clustering coefficient is considered. In our case, it is 
possible to obtain a more rich characteristic, namely, the 
distribution of the clustering coefficient in the graph. 

One can see that, in this graph, there is a one- 
to-one correspondence between clustering coefficient of 
a vertex and its degree: C = 2/k. Thus, the 
number m c (C,t) of vertices with clustering coeffi- 
cient C = 1, 2- 1 , 2~ 2 , . . . , 2 2 -*, 2 1 -*, 2-* is equal to 
3', 3 t ~ 1 , 3'~ 2 , . . . , 3 2 , 3, 3, respectively. In this case, it is 
natural to introduce the cumulative distribution of the 
clustering coefficient W cum (C) = J2c<c m c{C , t)/N t ~ 
c in3/in2 = (jy-l ^ w h er e C and C arepoints of the dis- 
crete spectrum. This corresponds to power-law behavior 



of the corresponding continuous distribution of clustering 
W(C) ~ C 7 ~ 2 for random scale-free network at small C. 

The average clustering coefficient can be easily ob- 
tained for arbitrary t, 
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For the infinite graph, C — 4/5, so the clustering is high. 

Degree correlations. — The number m(k,k',t) of edges, 
which connect vertices of degree k and k', characterizes 
short-range degree-degree correlations in the graph. It 
is convenient to write k = 2 P+1 and use the notation 
m(k,k',t) = c(p,p',t). Then one can find directly: 

c(t, t, t) = 3 , 

c(t,p' <t-l,t) = 3-2 t - 1 - p ' , 

c(p < t - l,p' < p - 1, t) = 3*~ p 2 p - p '- 1 . (3) 

This yields the cumulative distribution ~ fc 2 ~ 7 fc' _1 (we 
assume that k 3> k'), which, in turn, corresponds to the 
effective continuous distribution 



P(k,k') ~ k^k'' 



(4) 



This expression coincides with the corresponding asymp- 
totic formula for an arbitrary random scale-free citation 
graph |(| (by definition, a citation graph is a growing net- 
work, in which new edges do not emerge between pairs 
of old vertices) . Originally, Eq. ([|) was obtained exactly 
for a specific model in Ref. p2| . 

Shortest-path length distribution. — Here we briefly 
outline our exact results for the distribution V(£,t) = 
n(£, t)/[N t (N t — l)/2}, where n(£, t) is the number of pairs 
of vertices with minimal separation I. Details of the solu- 
tion and general expressions for n(£, t) will be published 
elsewhere. 

For calculation of n(£,t) one may use the following 
property. The length £^j of the shortest path between 
vertices i and j is equal to the minimal power of the ad- 
jacency matrix with nonzero {ij} element: {A 1 ^ 1 }^ — 
0, {A l }ij 7^ 0. This property allows us to obtain n(£,t) 
by counting the total numbers of nonzero elements in 
sequential powers of the adjacency matrix. This yields 
n(£,t) for several first steps: 
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9 6 

27 57 21 

81 351 369 60 

243 1806 3582 1716 156 



(5) 



and 



where t labels lines (t — 0,1,2,3,4,. 
1, 2, 3, 4, 5, . . . is the index of columns. 

The exact analytical form of the distribution V(£,t) 
was obtained by solution of recursion relations for 
n(£,t). In particular, an exact expression for the aver- 
age shortest-path length is of the form 
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£(t > 0) 



(4< + 11)3 
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10-3* 



3(3* + l)(3-3* + l) 



(6) 



One may check this expression using Eq. (^). The distri- 
bution quickly approaches an asymptotic regime, where 

^»l)-^+y + 0(t3- t ) = ^ln7V t + ^ln^. (7) 

Thus, the average shortest-path length logarithmically 
grows with increasing size of the graph. Expression 
(pp may be compared with the standard estimate p9[ : 
I ~ In N/ In k = In N/ In 4. The relative difference is sur- 
prising ly small, (4/9) / (In 4) - 1 = 0.616.... Notice that 
according to standard arguments the classical for- 

mula is not applicable for 7 < 3. Nevertheless, Eq. (Q) 
demonstrates that the classical estimate is unexpectedly 
good in this case where degree-degree correlations are 
strong. 

At large t, the distribution takes the Gaussian form 



P(£,t) = 



v /2 7 r(2 2 /3 3 )t 



exp 



2(2 2 /3 3 )t 



(8) 



which is violated only in narrow regions of width ~ i 1 / 3 
near the points £ = 1 and £ = t + 1 . One sees that the 
width of the distribution is of the order of \ft ~ \/ln Nt <C 
I(t). Notice that the simulations of the Barabasi- Albert 
growing random network also yield a Gaussian-like V{€) 



Eigenvalue spectrum of the adjacency matrix. — The 
spectrum G(A) contains N t eigenvalues. For t > 2, 
Nt-i—3 of them are equal to zero and, for t > 3, there are 
N t -2 — 3 eigenvalues equal to \/2 and the same number 
of those equal to — \/2. Here we do not derive analytical 
results for the entire eigenvalue spectrum but only study 
its tail using a simple numerical analysis. 

It is convenient to consider a cumulative distribution of 
eigenvalues G cum (A) = J2\>>\G(\'). The results of nu- 
merical diagonalization of the adjacency matrix for sev- 
eral timesteps are shown in Fig. [|. One sees that, in the 
large graph limit, the resulting cumulative distribution 
approaches a staircase-like form for A ^> 1. 

We found network-size-independent points of 
the spectrum. Using the coordinates of these 
points, we calculated the series of slops of 
lines connecting these points in Fig. |^: 
1.49847, 2.38192, 3.03023, 3.40683, 3.53135, 3.557 
actually a very good estimate ob- 
Interpolation of these values yields 

1 = 3.575 ± 0.015 of the cumula- 
^-(5-1). This value is 

. . to within a precision of 
our numerics. The height of steps in Fig. |^ is log 
their width approaches [(2/log 10 3) + (l/log 10 2)] 
A>1. 

Thus, one can suggest that exponent S of the corre- 
sponding continuous eigenvalue spectrum G(A 1) ~ 
A~ 5 is 8 = 2 + 7. One should mention that the direct 



(the last value is 
tained for t = 
the exponent <5 — 
tive distribution G cmrl (A) ~ 
2 + ln3/ln2 = 1+7= 3.585. 
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study of the eigenvalue spectrum for a growing random 
network with 7 = 3 (the Barabasi- Albert model) showed 
power-law dependence in a too narrow range of A to make 
precise conclusions |^3| (see also Ref. [Q). An estimate 
for the exponent in this situation was 6 ~ 5 p3| , which 
supports our conjecture. 

Percolation properties. — Let us delete at random ver- 
tices or edges of the pseudofractal. One may check that, 
in the infinite network limit, the giant connected compo- 
nent of the graph disappears only if "almost all" vertices 
or edges are removed. In the language of mathemat- 
ical graph theory this means that the fraction of ver- 
tices (edges), which we have to delete, approaches one 
as Nt — > 00. This is a standard property of scale- free 
networks with 7 < 3 [ [Tof . 

Discussion. — The network that we study in this Letter 
is a citation graph. Each new edge connects a new vertex 
and an old one. This is actually a deterministic variation 
of the scale-free growing network |25| ] in which one vertex 
is created per unit time and connects to both the ends of 
a randomly chosen edge. Therefore, it is not so strange 
that the properties of the pseudofractal network resem- 
bles those of scale-free random citation graphs. However, 
it is really surprising how close they appear to each other. 
Hardly one can propose more simple deterministic scale- 
free growing network. Therefore, we have found a very 
convenient tool for exploration of complex scale-free net- 
works. 

The extreme simplicity of the pseudofractal graph has 
allowed us to obtain a number of new results for grow- 
ing networks. In particular, for this network with strong 
correlations, we have obtained the exact shortest-path- 
lcngth distribution and the eigenvalue spectrum of the 
complex adjacency matrix. From the latter, we have 
made a conjecture that exponent of eigenvalue spectra 
of scale-free citation graphs is 2 + 7. Actually, we have 
shown that randomness of scale-free growing networks is 
of secondary importance for their structure. 
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FIG. 1. Scheme of the growth of the scale- free pseudofrac- 
tal graph. The growth starts from a single edge connecting 
two vertices at t = — 1. At each time step, every edge gen- 
erates an additional vertex, which is attached to both end 
vertices of the edge. Notice that the graph at timestep t + 1 
can be made by connecting together the three t-graphs. 




FIG. 2. Structure of the adjacency matrix of the graph 
(t = 2, Nt = 15). Black regions are unit elements of the 
matrix. In white regions, all matrix elements are zeros. In 
grey regions, nonzero (unit) elements are present. The ma- 
trix is symmetric, and each column in grey blocks above the 
diagonal contains only two nonzero elements. 
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FIG. 3. Log-log plot of the cumulative distribution of 
eigenvalues of the adjacency matrix, G cum (A) = J2y>\ C(A'). 
The curves show the spectra for t = 5, 6, 7, 8. The dashed 
line depicts the t — > oo limit. The t-independent points are 
marked. 
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